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Abstract 

Sivashinsky 's (1977) nonlinear integro-differential equation for the shape of corrugated 1- 
dimensional flames is ultimately reducible to a 2./V-body problem, involving the 2N complex poles 
of the flame slope. Thual, Frisch & Henon (1985) derived singular linear integral equations for the 
pole density in the limit of large steady wrinkles (N 3> 1), which they solved exactly for mono- 
coalesced periodic fronts of highest amplitude of wrinkling and approximately otherwise. Here we 
solve those analytically for isolated crests, next for monocoalesced then bicoalesced periodic flame 
patterns, whatever the (large-) amplitudes involved. We compare the analytically predicted pole 
densities and flame shapes to numerical results deduced from the pole-decomposition approach. 
Good agreement is obtained, even for moderately large Ns. The results are extended to give hints 
as to the dynamics of supplementary poles. Open problems are evoked. 
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I. INTRODUCTION 

Being able to describe the nonlinear development of the Landau-Darrieus (LD) 
instability of premixed-flame fronts is a central topic in combustion theory. As early as 1977 
Sivashinsky [3j showed, in the limit A <C 1 of small Attwood numbers based upon the fresh 
gas (p u ) or burnt gas (p& < p u ) densities, < A = (p u — pb)/ (p u + Pb) < 1, that the shape 
<f>(x, t) of a flat-on-average, spontaneously evolving wrinkled flame is governed by 

<k + \<fe = V<t>*x + I{<l>) (1) 

in suitable units. In ([1]) the subscripts denote partial derivatives with respect to time, 
t, and coordinate, x, normal to the mean direction of propagation, and the "viscosity" 
v > represents a reciprocal Peclet number based upon the actual flame thickness and 
the wrinkle wavelength. The linear integral operator /(•) is defined by I(e lkx ) = \k\ e lkx 
(whence I{4>) is the Hilbert transform, H(—(p x ), of — <p x ) and stems from the LD instability. 
The growth/decay rate of infinitesimal harmonics is \k\ — uk 2 , which identifies \jv and 
v as neutral wavenumber and minimum growth time, respectively. The nonlinearity is 
geometrical, accounting as it does for the cosine, (1 + s 2 ) -1 ^ 2 ~ 1 — s 2 /2 + . . ., of the small 
angle (arctan(s) ~ s + . . .) that the local normal to the flame front makes to the mean 
direction of propagation, where s ~ (f) x x A is the unsealed front slope. Originally derived 
m o as a leading order result for A —> + , equation ([!]) happens to govern the shape of 



steadily propagating fronts even when two more terms of the .A-expansion are retained 



Numerics [6] reveals that "steady" solutions of (CQ), corresponding to 4>(x,t) = —Vt + <j)(x) 
are often ultimately reached. When (CQ) is integrated with periodic boundary conditions for 
"not-too-small" values of v, v > 1/25 say, the "steady" pattern has a single crest per x- 
wise interval of 2ir length, where <p xx is large and negative; without loss of generality one 
may assume that one is located at x = 0, in which case <p x = when x is an integer 
multiple of tc (i.e., x = (mod 7r)) and <f) xx (±7r) ~ 1/tt. If Neumann conditions at x = 
and x = 7r are used instead, still with a moderately small z/, the final pattern obtained 
from numerical (pseudo-spectral) integrations of (0Q) may also have an extra crest located at 
x — 7T with (f) xx {Tr) large and negative. By the very way they are obtained as final state 
of an unsteady process the 2-crested patterns have a finite basin of attraction, contrary to 
the case of periodic boundary conditions [7] where the only stable patterns have a single 



its structure then remains valid practically up to A = 3/4, i.e., p u = 7pi, {4]. 
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crest per cell; yet such "half-channel" solutions happen to coincide with the restriction to 
< x < 7r of properly shifted 27r-periodic ones, for these are symmetric about x = and 
x = it. If v is too small the widest patterns get very sensitive to noise, even when caused 
by numerical rounding-off. In [8J the estimate fi > 0(e 1 / 2uK } = /i c (ft) was obtained for the 
noise intensity /i needed to trigger the appearance of extra-cells on top of the main ones with 
periodic boundary conditions; the number k in the above exponent is <f) xx (±7r) ~ I/71"; since 
;he most rapidly growing noise-induced disturbances (with initial wavenumbers \k\ ~ \ jv 
8j) of a nearly parabolic trough undergo an 0(e l l 2uK ) amplification, they ultimately get 
visible as subwrinkles of 0(1) final amplitude if /i > /U c (k). Having a larger (f) xx > 
at their troughs (see Sec. 7), 2-crested patterns are presumably less sensitive to noise 
than the single-crest ones associated with the same wavelength, because [i c (n) increases 
dramatically with k when v is small. The numerical work of Ref. 9fl also showed that sums 
4>(xi,x 2 ,t) = 4>i(xi,t) + (p2(%2,t) of orthogonal, 2-crested one- dimensional patterns play a 
central role in the study of (OQ) generalized to 2-dimensional flames (x — > (x\, x 2 ), 4> 2 X — > | V0| 2 , 
<p xx — > A0, /(•) = multiplication by [k ■ k) 1 ^ 2 in the 2-D Fourier space k = (fci,^)) 
and to rectangular domains in the Cartesian (xi,x 2 ) plane. Without noise such sums are 
exact stable solutions; with random additive forcing they recurrently appear as long-lived 
transients when Neumann conditions are adopted. 

Further analyses on the stability of solutions of Eq. ([1]) and their responses thus seem 
warranted, and getting the "steady" patterns that correspond to wide, hence large, cells (or 
small z/s) is a prerequisite. The present contribution is intended to do this. 

It is organised as follows. Section 2 introduces the pole-decomposition method, the 
discrete equations for the pole locations, and the two integral equations that approximate 
them for large front wrinkles. The latter equations are next solved analytically for isolated 
crests (Sec. 3) then one-crested periodic patterns (Section 4), and the prediction compared 
to numerical results from the pole-decomposition approach. Sections 5 and 6 compute the 
flame speed from the density, and take up the dynamics of a few extra-poles, respectively. 
Section 7 generalizes the above integral equations to a pair of coupled ones corresponding to 
2-crest periodic flames (and "half-channel" ones), then solves them analytically; comparisons 
with numerics are again presented. We end up with concluding remarks and open problems 
(Sec. 8). 
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II. POLE-DECOMPOSITION(S) 



In 1985 Thual, Frisch and Henon [lCj (herein referred to as "TFH") discovered (see also 



that ([I]) possesses solutions (p(x,t) representing 27r-periodic flame patterns with slopes 
(f) x in the form 

<f> x (x, t) = -uJ2 cot ( ) > ( 2 ) 

a=-N ^ ' 

in which the complex-valued poles of <p x (x,t), z a (t), are involved in conjugate pairs (z_ Q = 
z*, a ^ 0) for (f) x (x,t) to be real when x is. For this pole-decomposed expression to solve 
dl]), the z' a s (a = —N, . . . , —1, 1, . . . , N) must evolve according to the 2iV-body problem 

dz a N fzB-Z a \ 

-jf = v cot ( — ) ~~ 1 sign ( Jm ( 3 ) 

j3=-N ^ ' 

where Im (•) denotes the imaginary parts of (■) and the signum function (with sign(0) = 0) 
accounts for the LD instability. Once ([3]) is solved for the pole locations, 4>{x, t) is available 
from (T5]) and the wrinkling- induced excess propagation speed V = — (0 t ) > follows from 

V = I m , (4) 

where (•) stands for an average along the x-coordinate; thus, V simply measures the 
wrinkling-induced fractional increase in flame arclength, since ((1+s 2 ) 1 / 2 — 1) = (s 2 /2)+. . . ~ 



A 2 x V. Beside periodic 0(x, t)s, ([!]) also allows [10| for isolated non periodic wrinkles that 



have an infinite wavelength, V = 0, cot(z) replaced by 1/z, and 

dz a N 2 

— ^- = z/ V" isign(lm(z a )). (5) 

In the latter situation, the precise value of v > does not matter since it could be scaled 
out, and the integer iV > 1 is arbitrary. As for (J2]) Q^the maximum allowed value N opt (u) 
of N in steady configurations increases with 1/z/ > 1 [10]. As shown by TFH, steady flames 
obtained from (131) or ([5]) correspond to poles that "coalesce" (or align) along parallels to 
the imaginary axis, as a result of the pairwise pole interactions that are attractive along the 
real x-axis and repulsive in the normal direction. In the case of an isolated crest located at 
x = 0, the poles ultimately involved in steady solution are of the form iB a , —N < a < N, 
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a 7^ 0, with real B a s satisfying coupled discrete equations deduced from 

N 



V 



£ n W = sign(iU (6) 

^ D a — tip 



The authors of Ref. 



/3=-N 



lOj also evidenced that the larger the number N of pole-pairs in 



such "vertical" steady alignments, the smoother the involved poles are distributed along the 
B coordinate, with B a+ i — B a well smaller than B^. This suggested TFH to replace the 
discrete sum in (JB} (or its analogue deduced from Q) by an integral over the continuous 
variable B, with such a continuous measure that P(B) dB is the number of poles located 
between B and B + dB; a constructive definition of P(B) is specified in (|20|) . In this 
continuous approximation the steady versions of (J2j) ([3]) are amenable to singular Fredholm 
integral equations, specifically: 

dB' = sign(E) (7) 

in the non-periodic situations (an isolated wrinkle at x — 0), and 

uP{B') coth ( — - — J dB' = sign(5) (8) 

for the monocoalesced 27r-periodic cases (one single crest per cell, at x = (mod 2tt)). In 
(J7j) (jHJ) -B denotes the pole imaginary coordinate, and the Cauchy principal parts f ■ dB' 
stem from the condition (3 ^ a on the sums featured in ([3]) (jSJ). Consistent with their 
interpretation as pole densities, the P(B)s showing up in ([7]) (jHJ) both are non-negative even 
functions of their argument (for <p x to be real when x is) and are normalized by 

J P(B') dB' = 2N. (9) 

In ([7])- (ED the integrals extend over the ranges (to be determined as part of the solutions) 
where P(B) ^ 0. The next sections will solve (jTJ) (jSj) (jSJ) analytically, starting with the 
simpler equation (JTJ). 



III. ISOLATED CREST 

Because isolated crests have (fi x — ► at |x| — > cxo, we firstly anticipate the existence of 
some finite B mSuX > such that P{\B\ > B max ) = in (JTJ). We next recall the identity 



r /2 cos((2M+ l)<I) , )cos<I) , . 

T — ^ • L ^' = vrsm 2M+1 $ 10 



that can be deduced, through the change of variable $ — > $ + tt/2, from a similar one 
appearing in the Prandtl theory of lifting lines [12I . fl3 |. Identity ffTUl) allows one to solve 
such singular integral equations as Wigner's [ijj] (for the density, 2vP say, of eigenvalues of 
large real random matrices in the Gaussian Orthogonal Ensemble), written here as 

2vP(B'\ 



dB' = B- (11) 
B — B' y ' 

its solution is the celebrated semi-circle law 2txvP(B) = max(P max cos $, 0), [yj], provided 



that one sets 



B = 5 max sinf, -^<$<^, (12) 



in ffTTj) . Interestingly, the same change of independent variable in ([7]) produces 

"/ 2 2vP(&) cos $' . . , , 

^ — d$' = sign $ , 13 

. 7r/2 sin^-sin^' 5 v ; ' v ; 

since sign(P) = sign($) for |$| < tt. Over the same range (and hence over the narrower 
support of P, |$| < 7r/2), the right-hand side of ffl3|) may be expanded as the Fourier series 

.00 

sign($) = - 2M~J sin (( 2M + ( 14 ) 

^ M=0 

consistent with our convention that sign(0) = 0. From (flQl) the solution to (|T3|) can thus be 
written as a Fourier series of cosines that all vanish at $ = ±7r/2: 

4 00 1 

2 ^P($) = £ cos ((2M + 1)$) (15) 

= J_ log (I±^) (16) 

7T V 1 — COS $ / 



1 / 1 + Jl - B*/B, 



log v I ^ , (17) 



max 



and P = for |P| > B max ; to get ( 1171) from ( fl6l) . ( fl2j) was explicitly employed. 

The cumulative pole distribution -R(-B) = Jq P(B') dB' reads, after integration by parts, 

as 

2u R(B) = %^ f sin $ log + 2$) , (18) 

7T 2 V 1 - COS $ / 

whereby the renormalization condition R(B max ) = i?($ = tt/2) = iV fixes -B max to be given 
by 

P max = 2itNv. (19) 



TFH 10j fitted the cumulative distribution they obtained from a numerical resolu- 



tion of © for steady arrangements of aligned poles, by the expression it 2 vR 



f*\og(l.28Nv7r 2 /\B'\)dB' when \B\ < 5 max 10]. Equations flU} (PSD show that 1.28 
estimated from their numerical pole distribution at \B\ <C B max actually was a numerical 
approximation of 4/tt = 1.273 . . . Figures [T] and [2] compare the analytical findings (j!8p (IT9"|) 
to our own resolutions of ((SJ), with A" = 10, and 100, respectively. The TFH fit is also 
displayed for illustration. The pole density P is defined for a > 1 by 

P((B a + S a _i)/2) = {B a - B a ^)-\ (20) 

in terms of the pole locations (with B = by convention); it is shown in Fig. [3] for 
iV = 100, and compared with the continuous approximation (117j) and the TFH fit. Once the 
cumulative distribution is determined by (TLB"]) (1X9]) in the continuous limit, approximations 



B a to the discrete pole locations can be retrieved upon solving 



io| 



R(B a ) = a- 1/2, a = l,...,iV (21) 

numerically (e.g., by the Newton- Raphson method, with the "exact" B a s as initial guess!). 
The resulting crest shape 

N 

<f>(x) ■ 



-2^1og(l + £) (22) 



is compared to the exact one (numerical) in Fig. H] and to that issued from the continuous 
approximation. The latter profile has 

•B ma x 2uP(B) dB . , 

b (23) 

u max 



--sign(i) log ( ^===— | , (24) 



7T 



the second expression resulting from substitution of f fTTI) in ( l23l) . then a lucky look at p. 591 



of Ref.(l5] 



As suggested by the form of ( 1231) . and confirmed by (124"1) . x (x) is most simply deduced 
from P(±ix) through contour integration in the complex 5-plane. A further integration 
by parts of ff24|) yields the continuous-approximation prediction for 4>(x) (up to an additive 
constant): 

1 / cosh £ + 1 \ 

(j>{x) = -- sign(x) 5 max f sinh£ log cogh ^_ 1 + 2£j , (25) 



where x = -B max sinh£ (compare to (|T2l ). The integration constant was selected in Fig. 
H] to achieve good agreement with the exact 4>(x) for \x\ — > oo. Two final remarks: (i) v 
disappeared as a factor in (|2"11) as it should, because v can be scaled out; (ii) 0(x) is of 
the form vNF{x/vN), and this scale-invariance shows that the continuous approximation 
actually amounts to describing <p{x) at large distances compared to the actual radius of 
curvature (1/ J^ max AvP{B)dB / B 2 = o(is)) of the flame tip, when N ^> 1 ( that is, for large 
wrinkles) . 



IV. MONO-COALESCED, PERIODIC CREST 

The following simple remark will allow us to solve ([8]), i.e., in the case where all the 
poles of <p x are aligned along the imaginary x-axis (mod 2ir). Because P(B') still is an even 
function of B', only the even parts (at fixed B) of coth((5 — B')/2) will actually contribute 
to the integral over B'. Equation ([S]) may thus be re-written as 

f™ ' f ' 8 ''''-' rf ™ d B> = sign( B ), (26) 
J-B„„ tanh(B/2) - tanh(fl'/2) h 

upon use of the known formula for the tanh(-) of a difference, and neglect of a term propor- 
tional to f P(B') tanh(-B'/2) dB' = 0. We now set 

tanh(|)=tanh(^)sin$, -| < $ < |, (27) 

converting (Effil) into 

"" /2 2vP[B r ) cos $' . /JRN 

i — d® = sign($), (28 

n/2 sin^-sin^' B K h K ' 

which is nothing but ( 1131) . Therefore the sought after pole-density is still given by (1TB1) . 

the only difference with the previous non periodic case being that B, B max , and $ are now 

related by (J27D instead of (1121) . 

The new cumulative density = J P(S') dS' is given, after an integration by parts, 

by 

2 1, l + v4sin$ l + cos$ 

ti z vR{B) = -log log 

v ; 2 6 1 - A sin $ 6 1 - cos $ 



''" 1+Asin$' d$' 
1 - A sin $' sin $' 



+ / lo ST— T— ^7— S7, (29) 
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A = tanh(S max /2), whereby the normalization (jUJ) requires 

«r/2 l + Asin$ rf$ 
iVz/vr 2 = / log - 4 , ^ — — . 30 
7 1 - A sin $ sin $ v 7 

As the above integral turns out to be vrarcsinA (p. 591 of [l5() the range of P(B), still 
given by R(B max ) = N, now satisfies 

tanh( J B max /2) = sm(irNv) (31) 



instead of ( [19]) . The latter and ( !31j) coincide for ttNv <C 1, as do the associated pole densities. 
The maximum B max allowed by (13"Tj) . -B max = +oo, has 2iW = 1 and cos$ = 1/ cosh(_B/2), 
whence P(B) resumes the form 

P(B) = -L log (cotti^ (32) 



obtained by TFH wo Fourier transformations. The figure [5] compares our predictions (T29]) 
and (13 ip with very accurate solutions of ([3]) for iV = 100 and 2Nu = 1. Very good agree- 
ment is obtained even if N is only moderately large, and carries over to the pole densities 
themselves. Again, approximate solutions B a can be retrieved from the analogue of (|2ip . 
and an approximate flame front shape from 

N 

(j){x) = — 2z/^^log(l — cosx sech_B a ) + const. (33) 

a=l 

Figure shows of a comparison between (I33p . the exact flame shape obtained from the exact 
(yet obtained numerically) B a s satisfying and the curve deduced from the continuous 
approximation, for which the flame slope 4> x (x) reads 

J cot (^^) P(B) dB, (34) 

again a real function because P(—B) = P{B). With P(B) given by (TTB]) (l2"7j) (13"Tj) the above 
integral can be reduced to one available in p. 591 of 15| and yields (for — 7r < x < 7r): 

/ / \ 1 • , cosh£ + 1 x . . N 

0*0*0 = sign(f ) log — — — -, tan - = A smh f (35) 

7r cosh 4 — 1 2 

thereby confirming that <f) x (x) is accessible from P(B) by analytical continuation to ±ix. In 
particular, the TFH solution, eq. (I3"2j) . has = — 2 sign(x) log|cotx/4| and (^(±71") = 
1/tt; more generally, xx (±7r) = A/it . A further integration by parts yields 

. , , s . , N . 1 + sinh £ cosh £ + 1 

- 27T0(x) = sign(x) log — — — log — - 

1 — lA smh 4 cosh 4 — 1 

n ■ ^ 1 + iAsinhC dg 

*TZu&^*h- (36) 



-v 



which cannot be evaluated in simple closed form, but may be compared to (1291) ; of course 
4>(x) is real when x is, since the complex log(-) in fl36l) also reads 2iarctan(74sinh£) = ix. 
Note that 4>(x) has the form F(x; Nu), in the present units where the pattern is 27r-periodic. 
Adopting A ^ 2ir as wavelength would give 2n(j) = AF(2ttx/A] 2kNv/ 'A) with the same F. 
Accordingly, if uN/A is kept fixed, xx (±A/2) scales like 1/A as it should for v — > 0, whereby 
halving the wavelength renders the patterns less sensitive to noise (see the Introduction) . 



V. FLAME SPEED FROM CONTINUOUS POLE-DENSITY 



Plugging flMj) into (j4|) allows the wrinkling-induced increase in flame speed V to be 
written as 

2V = v 2 J J P{B) P{B') (cot cot x ~* B ) dB dB'. (37) 

Although the one-variable integrals involved when squaring ( 13^1) are ordinary ones, they 
may be written as principal parts. We next invoke the trigonometric identity cot a cot b = 
— 1 + cot (a — b) (cot b — cot a) and the average 

(cot ^—) = 1 sign(5) (38) 

to transform ( 1371) into 

2-^ = - J J P(B) P(B') dBdB' 
+2 jj sign(B)P(B) P(B') coth ( B ~ B j dB'dB (39) 

The first double integral (= (J P(B)dB) 2 ) in (|39l) follows from the normalization (Q, 
and is (2A) 2 . The second one is obtained from (jSj) after multiplication of both sides by 
P(B)sign(B)dB and subsequent integration over B: by ([9~lh it is 2N/v. Thus the simple 
formula 

y = 2Nv{l - Nv) (40) 



ensues; notice that it was obtained without having to solve ([8]). Actually (1391) can be shown 
from ([3]) to hold whatever N and z/ [161], again without solving the pole-equations themselves. 

In view of the accuracy of fj4*0|) one may inquire whether the solutions of ([7j) (jHJ) satisfy 
the "inviscid" Sivashinsky equation, i.e. (pQ) with v — 0, in the steady cases. To show they 
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do, for x 7^ at least, one may set V = Pv and M = Nv to remove v from ([7]) (Q, then 
process the Landau-Darrieus term of ([1]) as follows in the case of an isolated crest: 

/ x — iB 



2V(B) AJD [2V(B> 
x — iB 

4iV{B)V{B' 



f 2V{B') + < B ^ B <\ 
J B-B' y 1 

dB dB' = i<g, (41) 



(x — iB)(x — iB' 

where the notation (B' <-> B) represents a second copy of the integral that precedes it, with 
B and B' interchanged. The lines above successively use (I7|), acknowledge that (B,B') are 
dummy variables of integration that may be interchanged, then employ (1231) squared. Hence 
( |25|) satisfies ([1]) when v = and M is prescribed, if x ^ 0. Thanks to (1391) . a similar analysis 
applies to ([8]), provided i^O (mod 27r). 

Beside providing one with an exact P(B), eq. (Fl6l) shows that ([8]) admits a continuum 
of solutions, for there exists nothing in ([9]) to tell one that N ought to be an integer; this 
will be commented later (see Sec. 8). One finally specializes (JSj) to B = -£> max to show that 
./V is constrained by < 2Nu < 1, since coth(5 max — B) > 1 (see also (1311). 



VI. DYNAMICS OF SUPPLEMENTARY POLE-PAIRS 



In 2000, Vaynblatt & Matalon 13] addressed the linear stability of pole-decomposed 
monocoalesced "steady" solutions — Vt + (f>(x) to ([T|). Upon writing <fi(x,t) + Vt — <p(x) ~ 
exp(u;t)^„,(a:) <C 1 then analytically solving the linearised dynamics to get ui and ipu(x), the 



authors of 



17l | identified two types of linear modes. The modes of type I describe how the 2N 



poles of <fi x {x) evolve when displaced by infinitesimal amounts from equilibrium; all those are 
stable {oj < 0), but one that has to = (see below). The modes of type II were interpreted 
as resulting from x-periodic arrays of poles at ±ioo that may spontaneously approach 
the real axis if N is too small for the selected v < 1. The overall conclusion was thus: when 
endowed with 27r-periodic boundary conditions, all the monocoalesced solutions are linearly 
unstable, except a single one that has iV = A^ opt (z/) = [(1 + l/v)/2\ ~ l/2z/ (|_-J = integer 
part) and is neutrally stable (uj = 0) against shifts along the x-axis, the corresponding anti- 
symmetric eigenmode being ipo{ x ) = 4>x{%)- For iV < iV pt, modes of type II can manifest 
themselves, two particularly dangerous ones corresponding to incipient secondary wrinkles 
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centred on the main crests (x = 0, mod 2tt) or troughs (x = ir, mod 2tt). 

When Neumann conditions are employed instead, the aforementioned shifts are not al- 
lowed any longer because ipox 7^ at x = and x = ir. Numerical integrations 7j of ([T]) and 
(13]) evidence that there may then exist stable bi-coalesced patterns comprising an extra crest 
located at x = ir. Even though the steady 27r-periodic patterns also satisfy ([T]) with Neumann 
conditions when properly shifted to have 4> x (0) = = (j) x (n) no stability analysis similar to 



171 ] is yet available in this case; yet instabilities then necessarily require N < N opt (i>). Here 



we address a restricted aspect of the problem, namely: we study how the previously deter- 
mined monocoalesced "steady" solutions ( 1251) ( |36|) interact with extra pairs of poles. Since 
the free dynamics conserves the total number of pole pairs at its t = value, it makes 
sense to consider <f>(x, 0)s that involve them in a larger number (N + n) than the N = 0(1/ v) 
ones retained in a steady profile <p(x). Each of the n supplementary pairs at x m (t) ± iy m (t) 
contributes a perturbation (f> m (x,t) = (0 m ) — Au X^>i ex p( — j \Vm\) cos(j(x — x m ))/j to the 
flame shapes (this follows from (j2j) via a term-by-term Fourier expansion [10]) and, as shown 
in [181] , superposing m s can reproduce virtually any disturbance <p(x, 0) — 4>(x). In the 
present formulation the only difference between Neumann and 27r-periodic boundary condi- 
tions deals with the initial phases x m (0): whereas the former require the x m s to be compatible 
with the x <-> — x and it — x <-> it + x symmetries, the latter do not. 

Contrary to the more conventional normal-mode method (to which it is equivalent if 
|?/m(0)| 3> 1 ll8[), the pole approach can follow the disturbances when significant nonlinear 
effects set in. . . if one is able to solve the 2N + 2n coupled equations for the pole trajectories 
in the complex plane. The next remark somewhat simplifies the task. In the limits N ^> 1, 
v — > + and uN = 0(1) that led to ([8]), accounting for n = 0(1) extra pole pairs - as 
is assumed here - exerts only a small 0(v) perturbation on the 2N poles already aligned. 
Accordingly the distribution P(B) of poles along the main alignments at x = (mod27r) 
may be kept unchanged, and given by ( fTBT) ( |27l) ( |3TI) . when computing the motion of 2n 
supplementary ones. 

In the illustrative examples that follow only two extra poles (n = 1) located at 
±iy(t) (mod 2n), y > 0, then at ir ± iy(t) are considered, to begin with. 
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A. Extra-poles at x ~ (mod 2ir) 

When the two supplementary poles are located at ±iy(t) (mod 2ir), their altitude y(t) is 
determined from Q - within 0(v, 1/N) fractional errors in the limits iV ^> 1, v — > + and 
vN = 0(1) - by the ODE 

%=j vP(B') coth {^f-) dB> - 1, (42) 
2 

= — arcsin(sin(7riW) coth(y/2)) — 1, \y\ > B max , (43) 

7T 

where P(B) is the same as given by ffT6]) (|27|) and (l3~Tj) . to leading order, and leads to 
the closed form ( 1431) on integration [15]; for \y\ < -B max , dy/dt = by (jHJ). Therefore, 
whenever < 1 — 2^iV = 0(1) and v — > + , any initial y(0) > -B max will ultimately lead 
to y(+oo) = -Bmax; thereby adding one new incomer to the already present continuum. Put 
in words: if 2vN < 1 initially, the main pattern is unstable to disturbances with poles at 
±iy(t), and the latter process tends to make 2Nv approach 1 from below. 

Periodic boundary conditions would allow the supplementary pair to be initially off the 
x-axis, say at x(0) ± iy(0) with < x(0) < n (mod 2tt). The "horizontal" attraction by 
the main pole condensation at x = (0(u), actually) [lOj will make x(t) decrease, while 
y(t) still does if 2Nu < 1. Ultimately, the extra pole pair will join the main pole alignment 
(in finite time), and the previous conclusion is qualitatively unchanged: the process makes 
iV increase by one. When Neumann conditions are adopted, however, at least two pairs 
±x(t) ±iy(t) are needed if x(t) ^ 0, to meet the requirement of symmetry about x — 0, and 
two possibilities are encountered as to their fate. In the first instance, corresponding to not- 
too-small x(0)s and moderate values of y(0), the process is qualitatively the same as above, 
except that 2 pairs simultaneously join the main condensation at \y\ < B max , thereby making 
N increase by 2. If x(0) is small and y(0) well above B m3uX , the horizontal mutual attraction 
between the pair members may make them hit the x = axis at such a finite time t c that 
y(t c ) > -Bma X ; this is best shown from (jHJ) specialized to x(t) 1, whereby dx/dt ~ —vjx 
then x 2 (t) + 2i/(t — t c ) ~ for t < t c . The double pole thus formed at iy(t c ) then instantly 
splits into two simple ones lying on the x = axis at y(t) — y(t c ) ~ ±(t — tc) 1 ^ 2 , leading 
to a subsequent dynamics that ultimately ends like at the beginning of this subsection if 
2Nu < 1. 
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B. Extra-poles at x ~ ir (mod 2ir) 

In case the supplementary poles are located at tc ± iy(t) eq. (j4"2"l) is replaced by 

-j- = — arcsin(sin(7riW) tanh(?//2)) — 1 + z/cothy, (44) 
at ir 

since coth(w + Z7r/2) = tanhw. Even though z/ 1 the last term in (144)) . stemming from 
the interaction of the extra-pole with its complex conjugate, cannot be simply discarded, for 
otherwise (I44p would not be uniformly valid if y gets small. According to f)44p . any initial 
y(0) indeed ultimately leads to y(+oo) = v + o{y) and to a small (0(z/)) stable disturbance 
centred at x = ir (mod 2tt) whenever 2Nv < 1. Although the main pattern's curvature 
(pxxi 71 ) > is 0(1), and the 0(u) contribution to (p(x,t) of the extra pole-pair is small, it 



is nevertheless enough [7| to render the flame shape 4>{x,t) concave downward at x = ir; as 
shown in this occurs as soon as the extra poles enter a thin strip about the real axis, 
\y\ < V Aixu. Incorporating 0(N) extra pairs will also do, but the process of dynamical 
trough splitting is not within reach of such ODEs as (jHJ) when n = 0(N). The structure 
of 2-crested steady patterns with n = 0(N) will be studied in Sec JVIII 

Like in IVI Al one might begin generalizing the present discussion by envisaging a single 
pair of extra poles off the x = n axis, but this is already covered in the preceding paragraphs: 
if < a;(0) < 7r the pair ultimately joins the poles at x — ( mod 27r). It is more revealing to 
consider two such pairs at TT±x(t) ±iy(t) with x(t) "small enough" , in a way compatible with 
Neumann conditions, because something new appears. Comparatively large x(0)s will clearly 
lead to pairs that ultimately stick at x = (mod 27r), because their mutual horizontal 
attraction could not oppose that of the main alignments. The other extreme of very small 
a;(0)s again leads to the formation of double poles at some 7r ± iy(t c ), then a subsequent 
evolution of the two pairs 7r ± iyi^ along the line x = n (mod 2tt) until they settle at 
0{v) distances to the real axis if 2Nv < 1. The important conclusion is that stable 2-crest 
patterns exist when Neumann conditions are used and 2Nv < 1. 

By continuity there exist separating trajectories S±, such that none of the above be- 
haviours is observed if the pole pairs initially sit on them. The lines S± lead the two supple- 
mentary pairs towards an unstable equilibrium, a result of a competition between attraction 
by the main pole population at x = (mod 2ir), and the mutual attractions/repulsions 
among the pair members. For !/ < 1, and Nu = 0(1), using the steady version of ([3]) 
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and the pole-density given by ffT6l) (1271) (I3TT) . one can show that such equilibriums corre- 
spond to x(+oo) = ±(2nv / A) 1 / 2 + ■ ■ ■ and y(+oo) = ±u + ■ ■ ■ to leading order, again with 
A = tanh(i? max /2) = sin(7riW). This shows that there exist even more general steady so- 
lutions than considered elsewhere in the paper and in the literature (except in [7| where a 
similar conjecture was made on a numerical basis). One could have included other pairs as 
well, some of which along the x = tt (mod 2tt) axis. 

Our last remark is to again stress that the free dynamics ([3]) conserves the total number 
of poles (if finite). By the same token, allowing this number to vary with time is a means to 
study a forced version of the Sivashinsky equation: adding a pair of poles x m ± iy m aXt = t m 
amounts to accounting for a term 4> m (x) S(t—t m ) in the right-hand side of ([T]), and combining 
many (p m s with various phases (as to vary their signs), amplitudes (~ — 4z/exp(— \y m \) if 
\y m \ ^> 1) and times of implantation (t m ) could help one investigate the response of flames 
to a rich class of weak random noises. We understand that a similar proposal was developed 
about the "kicked" Burgers equation [19|], i.e., (DO) without the integral term in the one- 
dimensional case. 



VII. BI-COALESCED PERIODIC PATTERNS 



We now take up the structure of "steady" 27r-periodic solutions of (CD) that would have 
N pairs of poles i B a ( mod 2tt), a = ±1, ±2, . . . , ±iV, and n = O(N) other pairs at tt + i 6 7 
(mod 2tt), 7 = ±1, ±2, . . . , ±n. For brevity, we will say that the pole-alignments reside 
"at" x = or x = tt, respectively, like the two crests per-cell they correspond to. Because 
coth(w + z7r/2) = tanh(u), the steady versions of Q corresponding to such bi-coalesced 
flame patterns read as 




t ~* (^) 

5=—n ' 
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N fb — B \ 

+u ^ tanh ( 7 2 = si S n ( & 7)- ( 46 ) 



f3=-N 



In the distinguished limits v — > + , Nv = 0(1), nv = 0(1), the poles at x = and 
x = it get densely packed (at the scale of the wavelength), with densities P{B) and p(b), 
respectively. Both P and p will be nonnegative and, in general, compactly supported: 
P (\B\ > B max ) — — p(\b\ > & max ). The ranges -B max and 6 max are to be found as part of 
the solutions to the continuous versions of (145]) and fj46|) : 

y coth g J dS 7 

+ J vp(b') tanh f^y^-J = si S n ( s )> ( 47 ) 



y i/pt&O coth ^ 

+ J uP(B') tanh f^T") = si Mb), (48) 
that are valid for \B\ < -B max and |fe| < 6 max , respectively. To restore some symmetry we set 

tanh(5/2) = Asin$, A = tanh( J B max /2) < 1, (49) 
tanh(6/2) = asimp, a = tanh(6 max /2) < 1, (50) 

in ( l47l) ( l48l) . then acknowledge that both P(-) and p(-) are even functions, which allows one 
to suppress some odd parts of the integrands, viewed as functions of $' (or ip') at fixed $ 
(or ip). Some cumbersome algebra ultimately transforms (|47[) and (|48|) into 

2vP{&) cos $' 



sin $ — sin $ 



G?$' + sin $ 



/" 2vp(ip') cosy/ , , . . 

/ I 4 2 2 2 , 2 * d( P = S1 S n ( $ ) ( 51 ) 

J 1 — A^ar sin </?' sin <P 



/" 2vp{<p') cosy/ , . 

■f dip + Aa sm ip 

J sin <£> — sin ip' 

2vP(&) cos $' 



A 2 a 2 sin 2 $' sin 2 y> 



d& = sign(y>) 



(52) 
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where all the variables ($,$'), {<p,<p') are now taken in the common [— 7r/2,7r/2] range. 
One may thus adopt a common notation for them, (a, a') say, in both (15T1) and (152]) and 
subtract the results to eliminate the sign(-) functions in the right-hand sides. This produces 
an homogeneous equation for the difference P(-) — p(-), of which one obvious solution is 
P — p = 0. Hence the important result: if P = p is indeed a viable solution, then 

P{B) = j(a = arcsin(^^)), (53) 

p(b) = j(a = arcsin(^^)), (54) 

where J(o) is the same function for both. The even </(■) function itself is then found from 
(EH) or (J52D to satisfy 

" n/2 2u J (a') cos a' , , „ 
— da' + Aa 

/ 7r / 2 2uJ (<r') sin a cos o 1 , , , . 
. 2 . 2 da = sign(a). 55 
-tt/2 1 — A^cr sm a' sin a 

Further changing the independent variable to 9, with 

[1 + Aa) sin a 

1 + Aa sin 2 a 

fortunately converts the seemingly hopeless (1551) into a form equivalent to the already solved 
Eq. ([28"]) . playing the part that the former $ did there (most easily shown by starting 
from (jZHJ))- Accordingly, the solution to ( l55l) is available in terms of the already found pole- 
density pertaining to the isolated crests, then the monocoalesced ones: from ( TTBl) one indeed 
has 

r, TV \ 1 1 1 + COS^ . 

2i/ J a =— log, 7j; 57 

1 — COS 

with defined in ( 1551) . As said earlier, eqs. (1531) (1541) . -P(-B) is immediately retrieved upon 

setting sin a = tanh(£>/2) coth(j3 max /2) in (1561) (157|) : same operation to get p(b) from J(cx), 

upon setting sin a = tanh(6/2) coth(6 max /2) in (157|) (1561) . 

The first step to get -B max and 6 max again is to compute the cumulative pole-densities. 

For example R(B) = P(B') dB' is computed as follows from (157|) (1551) : 

l + cos#'\ 2/1 cos $'d$' 
cos 0' J 1 - A 2 sin 2 $' 



sint^ = — — - - - 2 _ , (56) 



2^^)=/ i°g(^^) r :::.z, . (68) 



1 + A sin $ , 1 + cos 9 

loe log 

s l-/4sin$ & l-cosfl 



+ 2 | S^ l08 T3T^ (59) 
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again with the understanding that $ (or $') is viewed as a function of 9 (or 9') via (1561) . and 
conversely; (1581) is obtained from the definition of R(B) upon setting tanh(5/2) = y4sin$, 
and (1591) results from an integration by parts. The cumulative density pertaining to p(b) is 
obtained in the same way from (|5"TI) fl5"6l) . now thanks to tanh(6/2) = a sin the result is like 
(I5"9"j) . except for the substitutions A — > a, $ — > if, B — > b, R{B) — > r(b). The normalisations 
R(B max ) = N and r(6 max ) = n thus impose the two conditions 



that may be compared to the former equation (130|) . and reduce to it when Aa = 0. Although 
we could not compute the above normalization integrals in closed forms, this can be done 
numerically without difficulty to get A and a as function of Nv and nv (or conversely). Note 
that iV > n is equivalent to A > a. N > n also implies that R(-) > r(-) when both are 
evaluated at the same argument, Fig. [HI 

Before closing this section, it remains to compare the above predictions to direct numerical 
resolutions of ( H5l) ( H6l) . by the Newton- Raphson method. This is done in Figs [7] and [HI 
Figure [7] will hopefully convince the reader that both P(B) and p(b) can be expressed in 
terms of the single function J(&) given by f[5~Tj) . 

Now that the pole densities are available, one may try to compute the corresponding 
increase in flame speed, V, from ( 1471) ( |48l) without solving them (like in Sec. 5), to produce 



this simple formula reduce to ( |40l) if n = 0, and could have been deduced directly from the 
discrete pole-equations, without solving them. The sum N + n plays the part N did for 
monocoalesced patterns and, as is shown upon specializing ( 1471) to B = 5 max , has to satisfy 
2(N + n)p < 1. 

As mentioned earlier, the flame slope <p x (x) pertaining to the continuous approxima- 
tion^) can be obtained directly from the corresponding pole-density(ies) via an analytical 
continuation from the real B (or b) axis to ±ix. Using the same procedure here gives, for 




(60) 



(61) 



V = 2v(N + n)(l-(N + n)v)- 



(62) 
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< x < tt: 



4> x = sign(x — x) log 



cosh if) + 1 



(63) 



7T 



cosh?/' — 1 



if) = - — ' - 

A — atari 2 x/2 



(1 + Aa) tan x/2 



(64) 



for bicoalesced flames, x being the point where sinh 2 if) 



oo in (164]) and, therefore, (p x (x) 



0: 



x = 2 arctan 



tanh B max /2 
tanh6 max /2 ' 



(65) 



At the flame front trough, k = <f> xx (x) = 2{A + a)/7r(l + Aa) > A/tt: the corresponding 
critical noise amplitude /z c (k) needed to trigger the appearance of subwrinkles markedly 
exceeds that pertaining to monocoalesced fronts (see Sec. 1). Two extra pole-pairs initially 
placed at the points ±x±iv (mod 2tt) would stay there in unstable equilibrium. There exist 
separating trajectories S± passing through them, which delineates the basins of attraction 
of the main pole condensations at x = or x = tt. Only the trajectories of initially remote 
extra-poles that are close enough to S± will enter the 0(y/v) strip adjacent to the B = 
axis where their direct influence on the main pattern becomes visible [7]. As seen from the 
real axis, the process then manifests itself as extra sharp sub-wrinkles seemingly "emitted" 
suddenly at x ~ ±x (mod 27r) before travelling to one of the main cusps where they even- 
tually join a main condensation. The iV/n-dependent shape of such separating trajectories 
thus controls the fate of "supplementary" poles of whatever origin, initial conditions or forc- 
ing; this will be exploited elsewhere, though one can already confirm that stable 2-crest 
steady patterns with n = O(N) ~ 1/u exist if Neumann conditions are employed. With 
27r-periodic conditions these are unstable even if 2u(N + n) = 1, as is seen by considering 
initial conditions where the 2n poles are slightly shifted to the left of x = tt (mod 2tt): both 
crests will ultimately merge. 

Comparisons with accurate numerical resolutions of the pole equations (145]) fj46l) are 
again good, Figj9l For v = 1/199.5, iV = 80, n = 20 they yield x = 2.053973 whereas our 
prediction (16"5T) . with A and a iteratively obtained from the normalization conditions (1601) 
(EE]), gives x = 2.053888. Like R(B) and r(6), 0(x) cannot be obtained in closed form, yet 
is readily accessible numerically. Also, if A = a, elementary trigonometry shows that the 
predicted flame slope (163]) resumes the result (135]) . up to a two-fold reduction in x- and 
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VIII. CONCLUDING REMARKS & OPEN PROBLEMS 



The above analyses may convey the feeling that the pole densities obtained so far have a 
family likeness, which is true because they were all deduced from the solution ([TBI) pertaining 
to isolated crests via adequate changes of independent variable. Whereas (IT5|) itself basically 
follows from standard Fourier analysis combined with a lucky re-summation of the series thus 
obtained (fT5l) . it would be interesting to understand why the changes of variable (1271) then 
(|56|) work so well. Admittedly the integral equation (I26I) bears some formal resemblance 
with ([7j), which guided us to propose the new variable fl27|) ; but that introduced in fl56|) 
looks more strange to us, and was actually discovered by trial-and-error after the 'resolvent' 
integral equation flo^l) is obtained. Yet fl5T)l) unlikely solely comes "out of the blue". In 
effect, one may note that fl56|) is equivalent to tanh(/3) = tanh(/3 max ) sin 9 if one defines 
tanh 2 (/3 max ) = tanh(_B max /2) tanh(6 max /2) and sets tanh(/3/2) = tanh(/3 max /2) sincx, which 
clearly mirrors what was employed to map the monocoalesced periodic case onto the isolated- 
crest problem. Hence ( l56l) rests on the celebrated composition law for hyperbolic tangents 
( r i , r 2) '■ T\ * r 2 = (ti + r 2 )/(l + TiT 2 ). It would be interesting to know whether the associated 
group properties give access to still more general solutions to the Sivashinsky equation 
([1]). That the scale- invariant signum function featured in ([7]) (151) is left unchanged by the 
successive changes of variables also is a key property that traces back to the presence of 
the Hilbert transform H(—(f) x ) in ([I]): in fine, it expresses that the complex velocity about 
the flame is a sectionally-analytic function in the complex x-plane, which is indeed a robust 
statement for it is little affected by conformal changes of variables that would leave the real 
axis globally invariant. 

Normalizing P(B) to 2A brings about the grouping Nv and, as long as the integral 
equations (J7J) of the continuous approximation are concerned, there is no reason why A 
should be an integer. Thus, © (jSJ) effectively admit a 1-parameter continuum of solutions. 
The situation is - in a sense, analogous to the Saffman- Taylor problem of viscous fingering 
and related ones (see 20j and the references therein): when surface effects (here curvature) 
are omitted, a continuum of steady patterns is found. The equation ([T]) for flames is peculiar, 
however, because one knows from the very beginning that only a discrete set of steady 
monocoalesced solutions exist, corresponding to As that are integers less than a well-defined 
^-dependent value, A opt (z/). The Sivashinsky equation (CQ) thus offers the opportunity to see 
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how the WKB approaches to finger-width selection developed 20| for the Saffman- Taylor 
problem, or kin, can be transposed to the present system to obtain a quantization condition 
on Nu; for here "inviscid" solutions are now available and one knows the answer in advance. 
This analysis likely is a key step to study flame response to noise, but has not yet been 
completed. Because WKB approaches essentially look for solutions of a linearised equation 
in the form exp(i j x k(x') dx'), where k(x) = 0(l/u) depends on the "inviscid" solution, 
it is seen that obtaining the latter to leading (0(1)) order in v is not enough. Hence the 
question: how to compute the leading {0{v)l) correction to the flame profiles obtained 
above? Obviously this would require to better understand the nature of the continuous 
approximation leading to the integral equations (JT])-© or (H7I) (14*81) for pole densities. In 
this context one may perhaps adopt the - rather unusual - view point that the exact pole 
equations (jSJ), once specialized to z a = iB a and steady patterns, constitute Gauss-like 
quadrature formulae to evaluate ([7]) numerically. How to define a "best" way of choosing 
the pivotal values, i.e. the B a s, naturally leads 2jJ to the notion of orthogonal polynomials 
associated with the Sivashinsky equation (JTJ). In the case of Wigner's equation (TTTjl the 
Hermite polynomials are invoked 141 ]. but we are not aware of such mathematical analyses 
about (HJ and © ©. 

Next, we recall that 2-crested patterns studied in Sec. VII also belong to a continuous 
family of solution profiles, now indexed by two independent parameters Nv, nv. Even if 
N + n is assumed to be given by the optimum value N op t(i/) ~ l/2z/, there still remains the 
question of how N/ n is selected in numerical resolutions of (CQ) with Neumann conditions at 
x = 0, 7T. The ratio N/n can undoubtedly be chosen by the initial flame shape <f>(pc, 0). In 
the case of forced propagations, the noise intensity (//) might well play also a role, for one 
can imagine situations where exp(— -k/2v) <C \i <C exp(— tc/Av): the noise is then intense 
enough to break monocoalesced patterns (see Introduction), yet too weak to noticeably affect 
the more curved 2-crested patterns with N = n. 

To tailor a global criterion as to compare the 2-crest patterns and their response to noise, 
the following remarks could be of some use. Let us collectively denote the B a s and 6 7 s as B 
and b, respectively. The unsteady versions of (]45p fj46l) - the pole equations for bicoalesced 
patterns - may be re-written as 

§ = -™ § = -™ < 66 > 
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in terms of U(B, b) = V(B) + v(b) + w(B, b), with 



B n — B 



(67) 



V(B) = v^\B a \-2v 2 lo g sinh 

a a,(3<a 

w(B, b) = -2u 2 log (cosh 67 j . ( 08 ) 

7,13 ^ 



and an expression similar to ( 1671) for v(b); t' = tjv is time scaled by the shortest growth time 
of small-scale wrinkles (see Introduction). Accordingly, when the right-hand sides of ( 1SU)) 
are supplemented with statistically identical, independent random (e.g., Gaussian) additive 
forcings, the joint probability density of (B, b) will tend to a quantity ~ exp(— U(B, b)//i 2 ), 
where /i -C 1 characterizes the noise intensities. Because U ~ 1 in the small-z/ limit (since 
P(B) and p(b) are 0(1 /V)) the above exponential is strongly peaked about the steady 
solutions. One can thus think of employing the iV/n-dependent scalar U(B,b), evaluated 
at steady state, as an objective means to discriminate the various 2-crested patterns in the 
presence of forcing. The task of evaluating U in the continuous approximation has not yet 
been completed. Neither is the analysis required to handle situations where the poles are 
slightly misaligned... yet still symmetric about x = and x = it for compatibility with 
Neumann boundary conditions. 

One must finally stress that the present analyses did not exhaust all the possibilities of 
"steady" solutions of ([1]), even with 27r as minimal periodicity. The interpolating solutions 
discovered in 3,0] constitute another class, comprising (possibly many-) extra poles, nearly 
evenly distributed [l?], Q] along sinuous curves at a distance from the real axis. In our 
opinion such unstable equilibriums are also worth analyzing in detail for v — > 0, as are those 
mentioned in Sec lVIl and generalizations of (pQ) itself 23]. 

As an end to a numerical work on (CQ), with noise included in the right hand side [7], one 
of us concluded that ". . . it is likely that new analytical studies of the Sivashinsky equation 
should be possible: even if the equation is now almost 30 years old, many things remain to 
be explained" . The words still hold true. 
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Fig. 1: Numerical vs analytical cumulative pole densities, for an isolated crest with 
1/v — 19.5, iV = 10. If exact, the theoretical curve (dot-dashed line, eq. (fT8"j) ) would pass 
through the middle of the risers of the numerical staircase (solid line, eq. (jS])). The dashed 
line is the TFH fit 23 

Fig. 2: Same as in Fig. [TJ for 1/v = 199.5, iV = 100. Only the upper hull (solid 
line) of the exact staircase is shown, for readability. The dashed line is the TFH fit 24 

Fig. 3: Numerical ( (TSUI) , solid line) vs analytical ( f|T71) . dot-dashed line) pole densi- 
ties P(B) for an isolated crest with 1/V = 199.5, N = 100. The dashed line is the TFH 
fit 25 

Fig. 4: Shapes of an isolated crest with 1/v = 19.5, iV = 10: continuous approxi- 
mation ((1251). dot-dashed line), exact (solid line), and smooth approximation from eq. ( |22l) 
(dotted) 26 

Fig. 5: Numerical (solid line) vs analytical (eq. ff29l . dot-dashed line) cumulative 
pole densities for a monocoalesced periodic crest, for 1/v = 199.5, and N = 100 
(= N opt (v)). Only the upper hull of the numerical staircase is shown 27 



Fig. 6: Shapes of a monocoalesced periodic flame with 1/v = 19.5, N = 10 (= N opt (v)): 
continuous approximation ( (136]) . dot-dashed curve) vs exact result (solid line) and smooth 
approximation ( (1331) . dotted) 28 



Fig. 7: Cumulative pole densities R(B) (upper curves) and r{b) for a bicoalesced 
periodic pattern with 1/v = 600.5, iV = 200, n = 100: the solid and the dotted lines are 
from eqs. (153]) fl54|) and (l5Tj) fl56l) : the dashed and the dot-dashed ones are the upper hulls 
of the exact staircases (see Fig. [T|). As (N + n) = N opt (v), _B 2 oo = oo 29 

Fig. 8: Theoretical pole densities P(B) (resp. p(b)) plotted as dot-dashed or dashed 



25 



lines vs 9, eq. (1571) . with sincr replaced by (ta.nhB/2)/A (resp. (tanh6/2)/a). The solid 
and the dotted lines are the numerical pole densities. All are for a bicoalesced periodic 
flame with N = 200, n = 100, \ jv = 600.5 30 

Fig. 9: Shapes of a bicoalesced flame with \jv = 199.5, iV = 80, n = 20: exact 
(solid line) vs continuous approximation (from integration of (!63j) . dot-dashed) 31 
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FIG. 1: Numerical vs analytical cumulative pole densities, for an isolated crest with 1/v = 19.5, 
N = 10. If exact, the theoretical curve (dot-dashed line, eq. (fTBj) ) would pass through the middle 
of the risers of the numerical staircase (solid line, eq. ©). The dashed line is the TFH fit. 
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FIG. 2: Same as in Fig. [IJ for 1/u = 199.5, N = 100. Only the upper hull (solid line) of the exact 
staircase is shown, for readability. The dashed line is the TFH fit. 
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— - theoretical pole density (formula 16) 
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pole vertical position 

FIG. 3: Numerical ( (j20|) . solid line) vs analytical ( (fT7|) . dot-dashed line) pole densities P(B) for 
an isolated crest with 1/u = 199.5, N = 100. The dashed line is the TFH fit. 
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FIG. 4: Shapes of an isolated crest with 1/v = 19.5, N = 10: continuous approximation ( (|25j) . 
dot-dashed line), exact (solid line), and smooth approximation from eq. (|22|) (dotted). 
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FIG. 5: Numerical (solid line) vs analytical (eq. (f29j) . dot-dashed line) cumulative pole densities 
for a monocoalesced periodic crest, for 1/v = 199.5, and N = 100 (= N opt (v)). Only the upper 
hull of the numerical staircase is shown. 
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x coordinate 

FIG. 6: Shapes of a mono coalesced periodic flame with 1/v = 19.5, N = 10 (= N op t(u)): continuous 
approximation ( ([36]) . dot-dashed curve) vs exact result (solid line) and smooth approximation ( ([331) . 
dotted). 
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FIG. 7: Cumulative pole densities R(B) (upper curves) and r(b) for a bicoalesced periodic pattern 
with 1/u = 600.5, N = 200, n = 100: the solid and the dotted lines are from eqs. (|53p (|54p and 
(fSTj) the dashed and the dot-dashed ones are the upper hulls of the exact staircases (see Fig. 
QJ. As (N + n) = N opt (iy), B 200 = oo. 
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FIG. 8: Theoretical pole densities P(B) (resp. p(b)) plotted as dot-dashed or dashed lines vs 6, 
eq. (|5Tj) . with sinu replaced by (tanhB/2)/A (resp. (tanh6/2)/a). The solid and the dotted lines 
are the numerical pole densities. All are for a bicoalesced periodic flame with N = 200, n = 100, 
\jv = 600.5. 
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FIG. 9: Shapes of a bicoalesced flame with 1/u = 199.5, N = 80, n = 20: exact (solid line) vs 
continuous approximation (from integration of (|63p . dot-dashed). 
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